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Many microorganisms find themselves immersed in fluids displaying non-Newtonian rheological 
properties such as viscoelasticity and shear-thinning viscosity. The effects of viscoelasticity on swim¬ 
ming at low Reynolds numbers have already received considerable attention, but much less is known 
about swimming in shear-thinning fluids. A general understanding of the fundamental question 
of how shear-thinning rheology influences swimming still remains elusive. To probe this question 
further, we study a spherical squirmer in a shear-thinning fluid using a combination of asymptotic 
analysis and numerical simulations. Shear-thinning rheology is found to affect a squirming swimmer 
in nontrivial and surprising ways; we predict and show instances of both faster and slower swimming 
depending on the surface actuation of the squirmer. We also illustrate that while a drag and thrust 
decomposition can provide insights into swimming in Newtonian fluids, extending this intuition to 
problems in complex media can prove problematic. 


I. INTRODUCTION 

Self-propulsion at small length scales is widely observed in biology, common examples include spermatozoa reaching 
the ovum during reproduction, microorganisms escaping predators and microbes foraging for food [U [5]. While 
swimming at low Reynolds numbers is well studied for Newtonian fluids [3], an understanding of the effects of complex 
(non-Newtonian) fluids on locomotion is still developing. Many biological fluids, such as blood or respiratory and 
cervical mucus, display complex rheological properties including viscoelasticity and shear-thinning viscosity HIS]- 
A viscoelastic fluid retains a memory of its flow history, whereas the viscosity of a shear-thinning fluid decreases 
with shear rate. While it is important to elucidate how non-Newtonian fluid rheology influences propulsion at low 
Reynolds numbers because microorganisms swim through biological fluids possessing these properties, an improved 
understanding may also guide the design of artificial micro-swimmers [B] and novel micro-systems [7] exploiting these 
nonlinear flnid properties. 

Recent research has begun to shed light on the effects of viscoelasticity (see the reviews by Sznitman and Arratia [8] 
and Elfring and Lauga i), but much less is known about swimming in shear-thinning fluids at low Reynolds numbers. 
Dasgupta et al. cni measured a decreased swimming speed of a waving sheet in a shear-thinning viscoelastic fluid 
relative to a Newtonian fluid. In contrast, an asymptotic study of a sheet driven by small-amplitude waves showed that 
the swimming speed of a waving sheet remains unchanged in an inelastic shear-thinning fluid compared to a Newtonian 
fluid m- A recent experiment by Gagnon et al. |12] on the locomotion of the nematode Caenorhabditis elegans has 
also suggested that shear-thinning viscosity does not modify the nematode’s beating kinematics or swimming speed. 
In addition, numerical studies [13 El examined a variety of two-dimensional swimmers and showed that faster or 
slower swimming in shear-thinning fluids can occur depending on the class of swimmer and its swimming gait. The 
results were understood in terms of the fluid viscosity distribution surrounding the thrust and drag elements of the 
swimmer. By estimating separately the propulsive thrust and drag force on the swimmer, Qiu et al. |B] obtained a 
scaling relation predicting the swimming velocity of a single-hinge swimmer (a micro-scallop), which is enabled to 
move at low Reynolds numbers by shear-thinning rheology. 
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The question that emerges from recent literature is when (and why) does a swimmer go faster or slower in a 
shear-thinning fluid m- To address this question we study a canonical idealized model swimmer, the squirmer, 
in a shear-thinning fluid described by the Carreau-Yasuda model using a combination of asymptotic analysis and 
numerical simulations.We predict and show instances of both faster and slower swimming depending on the surface 
actuation of the squirmer. We also explore separately the effects of shear-thinning on the propulsive thrust generated 
by the squirmer and the drag force it experiences, and demonstrate that extending these findings to swimming in 
non-Newtonian fluids can prove problematic. 


II. THEORETICAL FRAMEWORK 

The hydrodynamics of spherical bodies propelling themselves with surface distortions, otherwise known as squirmers, 
was first studied by Lighthill [16] and Blake m- We follow this approach and model a squirmer with prescribed 
time-independent tangential surface distortions. The resulting slip velocity around the squirmer is decomposed into 
a series of Legendre polynomials of the form ue(r = a,9) = (^)) where Vi{9) = —{2/l{l -f l))Pl {cos9) with 

Pi being the associated Legendre function of the first kind and 9 the polar angle measured with the axis of symmetry. 

The coefficients Bi are related to Stokes flow singularity solutions. In a Newtonian fluid, the Bi mode (a source 
dipole) is the only mode contributing to the swimming velocity, and the B 2 mode (the stresslet) is the slowest decaying 
spatial mode and thus dominates the far field velocity generated by squirmers. Therefore, often only the first two 
modes, Bi and B 2 , of the expansion are considered [TB1120j . The ratio of the two modes, a = B 2 /Bi^ characterises 
the type of swimmer in a Newtonian fluid: a > 0 describes a puller, which generates impetus from its front end 
{e.g. the alga Chlamydomonas), whereas a < 0 represents a pusher, which generates propulsion from its rear part 
{e.g. the bacterium Escherichia coli), and the a = 0 case corresponds to a neutral squirmer which induces a potential 
velocity field. In a Newtonian fluid, the swimming speed of a squirmer Un = 2Bij3 [121 dZ], which is independent of 
the fluid viscosity because drag and thrust change equally with viscosity. Any modes other than Bi only modify the 
surrounding flow structure but do not contribute to the swimming speed of a squirmer. This simple picture, however, 
does not apply to squirming in a shear-thinning fluid, as we discuss later, where all modes can potentially contribute 
to the swimming velocity, and adding any other modes to Bi will nontrivially affect the locomotion of the squirmer. 


A. Shear-thinning rheology: the Carreau-Yasuda model 

Shear-thinning fluids experience a loss in apparent viscosity with applied strain rates, a property that results from 
changes in the fluid microstructure. As the rate of strain exceeds the rate of structural relaxation, one observes 
microstructural ordering in the fluid m- Here, we capture the change in apparent viscosity due to this ordering using 
the Carreau-Yasuda model for generalised Newtonian fluids [22]. The variation of viscosity with applied strain rate 
is given by 


?7 = ’7oo + (?70-?7oo) [1 + A?|7p] = , (I) 

where 770 and 7700 are the zero- and infinite-shear rate viscosities respectively. The power law index n characterises 
the degree of shear-thinning {n < I) and the relaxation time At sets the crossover strain rate at which non-Newtonian 
behaviour starts becoming significant. The magnitude of the strain rate tensor is given by I 7 I = (n/2) ' , where 
n = jij'jij is the second-invariant of the tensor. As an example, measured values for human cervical mucus can be 
fitted by the Carreau-Yasuda model with values 770 = 145.7 Pas, 7700 = 0 Pas, At = 631.04 s, n = 0.27 [HIES]. 

We non-dimensionalise the flow quantities taking the first mode, Bi , of the surface actuation as the scale for velocity 
and radius, a, of the squirmer as the characteristic length scale. The strain rates are scaled with uj = Bi/a and the 
stresses by goco, such that the constitutive equation takes the dimensionless form 

r* = |/3+ (1-/3)[1 +Cubin'] (2) 

where t is the deviatoric stress tensor, and dimensionless quantities are denoted by stars (*). The Carreau number 
Cu = LoXt is the ratio of the characteristic strain rate, defined by the surface actuation w, to the crossover strain rate, 
defined by the fluid relaxation l/Aj. The viscosity ratio is given by /3 = riao/rjo G [0,1]. 

It is evident from that when the actuation rate w is much smaller or much larger than the fluid relaxation rate 
1/A(, i.e. when Cu —> 0 or Cu —> 00 , the shear-thinning fluid reduces to a Newtonian fluid of constant viscosity tjq 
(dimensionless viscosity 1) or 7700 (dimensionless viscosity j3) respectively. Recalling that for a given surface actuation 
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the swimming speed of a squirmer in the Newtonian regime is independent of the fluid viscosity, we therefore expect 
the swimming speed of a squirmer in a shear-thinning fluid to converge to its Newtonian value in the limits Cu —> 0 
or Cu —>■ oo. Non-monotonic variation of the swimming speed with Cu is expected for any swimmer with prescribed 
kinematics and has been observed by Montenegro-Johnson et al. [13j for some two-dimensional model swimmers. In 
this study we employ both asymptotic analysis and numerical simulations to investigate these effects of shear-thinning 
rheology on swimming at low Reynolds numbers. 


B. Asymptotic analysis 

The deviatoric stress tensor r* in (§ is a non-linear function of the strain rate tensor 7 *. Assuming only a weak 
nonlinearity, we may uncouple the Newtonian and non-Newtonian contributions writing 

r*=7*+eA*, (3) 

with £ 1 as a dimensionless measure of the deviation from the Newtonian case (£ = 0). 

We observe that in the limits Cu = 0 or /3 = 1 , (H reduces to a Newtonian constitutive equation. Thus one 
may expect weakly nonlinear behaviour when the fluid relaxation rate is much faster than the surface actuation rate 
(e = Cu^ <C 1 ), or when the zero-shear-rate viscosity is very close to the infinite-shear-rate viscosity (e = 1 — /3 <C 1). 
Henceforth we shall work in dimensionless quantities and therefore drop the stars (*) for convenience. 


1. Expansion in Carreau number 

Expanding all fields in regular perturbation series in e = Cu^, we obtain, order by order, the constitutive equations 

To = 7o> 

Ti =7i + ^’^^Y^(l-/3)l7op7o, (4) 

hence A = (1 — /3) |7op7o to leading order in Note that the first correction to the Newtonian behaviour is 

linear in n, which points to a linear dependence of the swimming speed in n upon using (|^, elucidating the trend 
suggested by the two-dimensional numerical findings in Montenegro-Johnson et al. [14]. We also remark that this 
expansion is valid only when I 7 P is o(l) and is therefore not uniformly valid across all values of strain rates. 


2. Expansion in viscosity ratio 


Expanding in perturbation series with £ = 1 — /3 gives us, order by order, the constitutive equations. 

To = 7o> (5) 

Ti =7i + |-l+(l + Cu2|7on^|7o: (6) 

where in this limit A = | —1 + (l + Cu^lyoP) ^ | 7 q to leading order in ([^. Note that this asymptotic expansion 
is uniformly valid for all strain rates or Carreau numbers, which permits a full-range study of the non-monotonic 
swimming behaviour. 


C. The reciprocal theorem 

Stone and Samuel |24] demonstrated the use of the Lorenz reciprocal theorem in low-Reynolds-number hydrodynam¬ 
ics [25] to obtain the swimming velocity of a squirmer for a given prescribed surface actuation u'^ without calculation 
of the unknown flow field, provided one can solve the resistance/mobility problem for the swimmer shape (with surface 
S). Lauga [MlEl] then developed integral theorems extending this method for use with complex fluids. We use these 
methods in the subsequent calculations to obtain the swimming velocity of a squirmer in a shear-thinning fluid, which 
closely follow the formulation in Elfring and Lauga j^. 
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We represent the velocity field and the associated total stress tensor for a force- and torque-free swimmer with u and 
(T respectively. We consider the corresponding resistance problem in a Newtonian fluid to simplify the calculation of the 
swimming velocity. The resistance problem (denoted with a ‘hat’) involves the rigid-body motion with translational 
velocity U and rotational velocity ft, and the corresponding velocity field and associated stress tensor are represented 
by u and & respectively. Due to the linearity of the Stokes equation, we may write u = C- U, a' = 'T-U and 
T = —'R- ■ lA. Here, for compactness, both the translational and rotational components of velocity are contained in U 
and, similarly, the corresponding matrices contain both the translational and rotational terms. In weakly nonlinear 
complex fluids, the six-dimensional translational and rotational velocity of the swimmer lA is then given by 


u = 



{n-f)dS-e [ A ; V£ dH 
Jv 


(7) 


The integral over the volume of fluid V external to S in the equation measures the change in swimming dynamics 
due to the non-Newtonian behaviour of the fluid. For a spherical squirmer with axisymmetrical tangential surface 
distortions, there is no rotational motion and the translational velocity is given simply by 


U = / u'^dS’ - / A: (l + VG dF, (8) 

dTTJs 8 tt Jy \ 6 J ^ ’ 

where G = y (l-f is the Oseen tensor (or Stokeslet). The first term on the right-hand side is the result of 
swimming in a Newtonian fluid |24j . and the last term in equation contains the weakly nonlinear effect, which can be 
evaluated analytically in some special cases and can be computed in general by numerical quadrature readily. 


D. Numerical solution 

The numerical simulations of the momentum equations at zero Reynolds number with the Carreau-Yasuda consti¬ 
tutive relation Q are implemented in the finite element method software COMSOL. We use a square computational 
domain of size 500a x 500a, discretized by about 30000-50000 Taylor-Hood {P2 — PI) triangular elements. The mesh 
is refined near the squirmer in order to properly capture the spatial variation of the viscosity. Since slowly decaying 
flow fields are expected at low Reynolds numbers, a large domain size is important to guarantee accuracy. The sim¬ 
ulations are performed in a reference frame moving with the swimmer and the far-held (inlet) velocity is varied to 
obtain a computed zero force on the squirmer. In addition to comparing with the asymptotic analysis in this paper, 
we have validated our implementation against the analytical results of a three-dimensional squirmer in a Newtonian 
huid [161117] and a two-dimensional counterpart in a shear-thinning huid [14] . 


III. RESULTS AND DISCUSSION 

As a Hrst step we investigate the effect of shear-thinning rheology upon swimming speed by considering the small 
Cu regime and use ([^ to derive an analytical formula for the leading-order swimming speed ?7 of a two-mode squirmer 
(with Bi and i ?2 modes) 


^ = l + Cu^{l-f3) [1 + C^a^] , (9) 

where Ci = 0.33 and C 2 = 2.25 are numerical constants, and Un is the Newtonian swimming speed. In a shear¬ 
thinning huid we have n < 1 and /3 < 1 , and hence we hnd that this two-mode squirmer can only swim slower than in 
a Newtonian huid {U/Un < 1) in the small Cu regime. The two-dimensional numerical simulations in Montenegro- 
Johnson et al. m reported that a neutral squirmer [a = 0) swims slower in a shear-thinning huid, which is consistent 
with our analytical results for a three-dimensional squirmer in the small Cu regime, but we also hnd that the 
same conclusion of a decreased swimming speed holds for pushers {a < 0) and pullers (a > 0) as well. In contrast 
to the case of swimming in a viscoelastic huid, where the pusher and puller attain diherent velocities given the same 
magnitude of a m, © reveals that in a shear-thinning huid a pusher and puller have the same swimming velocity 
because the function for swimming speed is even in a; this asymptotic result is verihed by numerical simulations to 
hold for diherent ranges of Cu and /3 (as shown by the overlapping of the upper and lower triangles in hgure[^. 

To further characterise the variation of swimming speed, over the full range of Cu, we consider the asymptotic limit 
e = 1 — /3 <C I, aided by numerical simulations for larger values of e. Biological huids often have a small viscosity ratio 
/3 and hence e = I — /? is typically close to 1. In hgure la, we present the numerical results for a neutral squirmer. 
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FIG. 1. In (a) we show results from numerical simulations for a neutral squirmer (q = 0, •), puller (a = 5, A) and pusher 
(a = —5, T) for e = 0.99 and n — 0.25. (b) The non-monotonic variation in velocity is well captured by asymptotics for a 
neutral squirmer (solid line), and a pusher/puller (dashed line) at e = 0.1 and n = 0.25. 


pusher, and puller in the biological limit using the values e = 0.99 and n = 0.25 to emulate human cervical mucus 
[mils]. We demonstrate in the upper inset (neutral squirmer) and lower inset (pusher and puller) in figure la that 
the numerical solutions for the swimming speed ratio converge to the asymptotic solutions (solid line in the upper 
inset; dashed line in the lower inset) when £ —>■ 0. In hgure lb, the results are presented at £ = 0.1 and we note 
that all qualitative features of the impact of a shear-thinning fluid in the biological limit (£ « 1, figure la) on the 
swimming speed are well captured by the asymptotic analysis (when £ <C 1, figure lb) and as expected, the numerical 
simulations (symbols) agree very well with the asymptotic theory (lines) when e is small (figure lb). 

The non-monontonic variation of the swimming speed with Cu may be expected based on the asymptotic behaviour 
of the constitutive relation discussed at the end of section |II A[ To understand the variation more quantitatively, 
recall the form of A in © and observe from the integral expression for swimming velocity (|^ that at low strain 
rates, the non-Newtonian contribution A ~ ^Cu^(n — l)|7op7o vanishes as Cu —)• 0. At high strain rates, A ^ 
— 7 g -I- (Crtlyol)" ^ 7o; the hrst term, —7o) vanishes under the integration in ([^ [5], and the remaining term gives a 
non-Newtonian contribution that vanishes as Cu —>■ oo because for a shear-thinning fluid n < 1. The swimming speed 
therefore displays a non-monotonic variation with Cu, and since the speed decreases when Cu is small as shown by 
©, a minimum swimming speed may be expected to occur at intermediate values of Cu (the ‘power-law’ regime of 
the model), where the non-Newtonian effect is most significant. However, for a given swimming gait, if the actuation 
rate of the swimmer is small enough or large enough, the shear thinning fluid may appear to have no effect at all on 
the swimming speed. 

To understand the reduction in swimming speed, inspired by the qualitative descriptions given in Montenegro- 
Johnson et al. m, we look into the thrust and the drag of the swimming problem separately. 


A. Drag and thrust 


We separate the swimming problem into a drag problem (a sphere undergoing rigid body translation U inducing 
hydrodynamic drag), and a thrust problem (a sphere held fixed undergoing only tangential surface distortions thereby 
generating thrust). The superposition of these two sub-problems gives the entire swimming problem in a Newtonian 
fluid in the Stokes regime; this is obviously not the case in a shear-thinning fluid due to its nonlinear constitutive 
equation. However, by looking at the thrust and the drag problems, one may gain insight into the more complex 
non-Newtonian swimming problem HE]. 

We derive the expressions for drag and thrust in a shear-thinning fluid again via the reciprocal theorem approach 
(section H CI by utilizing the solution to the resistance problem in a Newtonian fluid. The drag force on a sphere 
moving with a velocity U in weakly shear-thinning fluid is given by Fd = — 67rU — |£ Jy A^ : (l -|- VG, where 
A^) is formed by the solution to the Newtonian drag problem. Similarly, the thrust force generated by a sphere held 


stationary with surface actuation u'® in a weakly shear-thinning fluid is given by Ft = /. 


, u‘ 


'dS-leJyA, 


(l -I- gV^) VG, where At is formed by the solution to the Newtonian thrust problem. 

One could expect a drag reduction when a rigid sphere is pulled with a constant velocity through a shear-thinning 
fluid since the fluid viscosity is reduced by the fluid straining motion. However, it is interesting to see in figurethat 
the thrust reduction caused by the shear-thinning rheology is larger than the drag reduction for a large range of Cu. 
This more severe reduction in thrust than drag then suggests slower swimming speeds compared with the Newtonian 
case, which correctly predicts the trend found by detailed calculations (figure]^. In addition, for very small or large 
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FIG. 2. In (a), the symbol A denotes the difference between drag, or thrust, in a shear-thinning fluid and a Newtonian fluid: 
the dash-dot curve represents the drag reduction in a shear-thinning fluid compared with a Newtonian fluid, while the solid 
and dashed curves represent loss in thrust for squirmers with a = 0 and a = ±5 respectively (n = 0.25, £ = 0.1). In (b) we 
show that the difference between drag and thrust is positive both for a neutral squirmer (solid) and a pusher/puller (dashed) 
in a shear-thinning fluid. All quantities are dimensionless 


values of Cu, the difference between the magnitudes of drag and thrust {Fd — Ft) vanishes as shown in figure [^, 
respecting the limits where the swimming speed should recover the Newtonian value (figure [^. 

Although conceptually intuitive, the drag and thrust decomposition is not complete as it neglects the contribution 
of non-linear products in the non-Newtonian stress to the full swimming problem, namely A yf At + Ajj owing to 
the non-linearity in the constitutive equation. We will give a counter-example in section |IIIB| below to illustrate a 
scenario when these intuitive arguments fail. 


B. Addition of other squirming modes 


The results from the detailed asymptotic and numerical analysis as well as the intuitive model for a two-mode 
squirmer seem to suggest that the shear-thinning rheology acts to hinder the locomotion. This raises the simple 
question of whether this conclusion still holds if other modes of surface actuation are present. The picture is clear for 
a Newtonian fluid: only the Bi mode contributes to swimming and adding other modes does not alter the swimming 
speed. However, can the shear-thinning rheology render other modes, typically not considered in the Newtonian 
analysis, effective for propulsion? We address these questions using the asymptotic and numerical tools developed in 
the previous sections. 

We first note that the mode alone leads to locomotion in a shear-thinning fluid in stark contrast to a Newtonian 
fluid as only the Bi mode has a non-zero surface average (see §)• Indeed any odd mode alone may lead to locomotion 
in a shear-thinning fluid (even modes alone do not swim by symmetry). We also find quite distinctive behaviour when 
the B^ mode is combined with other modes. In the Cu 1 regime, we can derive an analytical expression allowing 
us to predict the values of a and = B^/Bi for faster or slower swimming. To quadratic order in Cu, we And 


Yj— — 1 + Cu^ (1 — /?) -— ^ —^Ci [i -I- C2 (1 + C3C) o? + C 4 (CsC^ + CgC ~ 1 ) C] I 


( 10 ) 


where the additional numerical constants are given by C 3 = 
the swimming speed is even in a and we recover when 
swimming (U/Un > 1) will occur if 


0.51, C 4 = 0.70, C 5 = 0.18, and Ce = 1.66. Again 
: 0 as expected. From (10) we can predict that faster 


d^U 

dCu^ 


(a^C) 


> 0 , 

Cu=0 


( 11 ) 


in other words when the term in the square brackets in (10) is negative. In figure]^ we plot the level set curve below 
which faster swimming occurs in the small Cu regime; we find that this can only occur when ^ is negative for any a. 
For example, when a = 0 we must have C < —10.11, while a = ±5, C < —2.22 leads to faster swimming. 

In figure we show the variation of swimming speed for two swimmers with a and C chosen below the level set 
curve (the upper solid and dashed lines) and two swimmers with a and C chosen above the level set curve (the lower 
solid and dashed lines). We note that the swimming speeds of the faster swimmers in the small Cu regime experience 
a subsequent fall below the Newtonian value and then a rise above it as Cu increases, before asymptoting to the 
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FIG. 3. In (a) we show the level set curve below which faster swimming occurs for specific values of a and C, = B^jBi when 
Cu ^1. In (b) we show the swimming speed of a neutral squirmer (solid line) and a pusher/puller (dashed line) at values of ( 
chosen below the curve in (a) so that the swimming speed is larger than the Newtonian value at small Cu (the upper solid line: 
q: = 0, (( = —15 ; the upper dashed line: a = ±5, ^ = —4 ); conversely, with values of ^ above the curve in (a) the swimmers 
always swim slower than in a Newtonian fluid as shown (lower solid line: a = 0, C = 15; lower dashed lines: a = ±5, C = 4). 
Here e = 0.1, n = 0.25. 


Newtonian swimming speed at high Cu. This indicates that microorganisms (with a given swimming gait) can swim 
both faster, and slower, than in a Newtonian fluid depending on the actuation rate of that gait. In contrast, the 
two swimmers that swim slower in a Newtonian fluid in the small Cu regime remain slower for larger Cu with a 
non-monotonic variation similar to that observed previously (see figure [^. These results also hold qualitatively for 
large values of e. 

We emphasize that the thrust and drag reduction model is unable to explain the faster swimming speed with the 
addition of a mode because in these cases the thrust still decreases more than the drag over a wide range of 
Carreau numbers, if they are considered separately. This serves as a counter-example demonstrating how analyzing 
drag and thrust separately may not adequately describe swimming in complex media. This fact points specifically to 
the interaction between thrust and drag fields, due to the nonlinearity in the constitutive equation, as the cause of 
faster than Newtonian swimming. 


IV. CONCLUSION 

We show in this paper that shear-thinning rheology affects a squirmer with a prescribed swimming gait in nontrivial 
and surprising ways; we predict, analytically, instances of both faster and slower swimming than in a Newtonian fluid 
depending on the details of the prescribed boundary conditions. Indeed we demonstrate that even with the same 
squirming modes a squirmer can swim faster or slower depending on its rate of actuation. In general, these results point 
to the importance of both the spatial and temporal details of the swimming gait of a microorganism and ultimately the 
difficulty in predicting the resulting effect of the non-Newtonian fluid a priori. In light of this, an important next step 
would be to incorporate models of internal force generation for biological swimmers and determine how fluid rheology 
affects the resultant gait itself in concert with propulsion. Finally, we remark that the drag and thrust decomposition 
of the swimming problem is indeed effective in Newtonian fluids and may also be insightful in complex fluids in some 
instances, but one should use caution when extending the results to non-Newtonian swimming as the inherent non¬ 
linearity of the problem can be significant enough for a Newtonian-like decomposition to yield qualitatively flawed 
predictions as illustrated by the example we provide. 
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